On the Integrability of the Rabi Model 
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The exact spectrum of the Rabi hamiltonian is analytically found for arbitrary coupling strength 
and detuning. I present a criterion for integrability of quantum systems containing discrete degrees 
of freedom which shows that in this case a finite symmetry group may entail integrability, even 
without the presence of conserved charges beyond the hamiltonian itself. Moreover, I introduce and 
solve a natural generahzation of the Rabi model which has no symmetries and is therefore probably 
the smallest non-integrable physical system. 
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The Rabi (or single-mode spin-boson) model consti- 
tutes probably the simplest physical system beyond the 
harmonic oscillator. Introduced over 70 years ago 1], 
its applications range from quantum optics Q and mag- 
netic resonance to solid state [sl and molecular physics 
[3] • Very recently, it has gained a prominent role in novel 
fields of research such as cavity QED 5] and circuit QED 
0. It can be experimentally realized in Josephson junc- 
tions or using tr app ed ions Q, in Cooper-pair boxes 
and flux q-bits [lO]. In this way, its complete theo- 
retical understanding is mandatory for all feasible ap- 
proaches to quantum computing [ll|. Despite its old 
age and central importance, the Rabi model has not 
been exactly solved JJj, il3j] . It shares with the other 
paradigma of quantum physics, the hydrogen atom, an 
infinite-dimensional state space but — in contrast to the 
latter — the spectrum and eigenfunctions of the Rabi 
model are known only by numerical diagonalization in a 
truncated, flnite-dimensional Hilbert space. This is quite 
surprising, as the Rabi model has a smaller number of 
degrees of freedom than the hydrogen atom. In particu- 
lar, a single degree of freedom a;, subject to a harmonic 
potential, couples to a quantum system with only two 
allowed states | t) and | J,). Therefore, it does not pos- 
sess a classical limit: the quantum degree of freedom has 
a finite-dimensional Hilbert space and places the Rabi 
model in between the case of one and two (classical) de- 
grees of freedom. The hamiltonian reads in units with 



(1) 



Here, the Ux^z are Pauli matrices, describing the two-level 
system with level-splitting 2 A and a {a^) denote destruc- 
tion (creation) operators of a single bosonic mode with 
frequency These two systems are coupled through a 
term proportional to which has different interpreta- 
tions according to the experimental situation to model. 



Although ([T]) describes the simplest of all physically sen- 
sible interacting quantum systems, it poses a serious ob- 
struction to its analytical solution because of the appar- 
ent lack of a second conserved quantity besides the en- 
ergy, which has led to the widespread opinion that it 
cannot be integrable [l3 - [T9| . To remedy this difficulty, 
Jaynes and Cummings proposed already in the sixties an 
approximation to ([l]) which does possess such a quantity 
1201. Their hamiltonian reads 



AcT^ 



(2) 



fltfl -f \{<Jz + 1) commutes 
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Hjc = Lja'a + g{(T 

with (7^ = {ax ± i(Ty)/2. 
Here, the operator C 
with Hjc and leads at once to the solvability of ([2]). 
The Jaynes-Cummings model is the so-called "rotating- 
wave" approximation to ([T]) and was justified because 
the conditions of near-resonance 2A « uj and weak cou- 
pling g <^ Lo for such an approximation are realized 
in many experiments. The conservation of C means 
that the state space decomposes into an infinite sum of 
two-dimensional invariant subspaces, each labeled by the 
value of C = 0, 1, 2, . . . Each eigenstate of ([2]) is then char- 
acterized by two quantum numbers, namely C and a two- 
valued index, for example -I- and — , denoting a basis vec- 
tor in the two-dimensional subspace which belongs to C. 
Whereas the possible values of C form an unbounded set, 
corresponding to the quantization of a classical degree of 
freedom, the second quantum number can take only two 
values, refiecting the intrinsic quantum nature of the two- 
level system. The conserved quantity C generates a con- 
tinuous J7(l)-symmetry of the Jaynes-Cummings model 
([2]) which is broken down to Z2 in the Rabi model ([T]) due 
to the presence of the term a^cr^ -I- aa~ . This residual 
Z2-symmetry, usually called parity, leads to a decomposi- 
tion of the state space into just two subspaces each 
with infinite dimension. One would conclude that this 
symmetry cannot suffice to solve the model exactly — 
but in fact it does. Whereas a discrete symmetry is too 
weak to accommodate a classical (continuous) degree of 
freedom, it can do so with a quantum degree of freedom 
— like the two-level system in the Rabi model. 

We observe a direct relation between the nature of the 
degree of freedom (continuous or discrete) and the sym- 
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metry (C/(l) versus Z2), which can be used to "eliminate" 
it by fixing the corresponding irreducible representation. 

Our main result is the following: The spectrum of ([T]) 
consists of two parts, the regular and the exceptional spec- 
trum. Almost all eigenvalues are regular and given by the 
zeroes of the transcendental function G± (x) in the vari- 
able a;, which is defined through its power series in the 
coupling g: 



G±{x) = J2Kn{x) 



It 



A 



(3) 



The coefficients Kn{x) are defined recursively, 

nKn = fn-l{x)Kn-l - (4) 

with the initial condition Kq ~ l^Ki{x) — fo{x) and 

A2 



Jn{x) = h — [nuj- X 

UJ Zg \ X — nu! 



(5) 



The function G±{x) is not analytic in x but has sim- 
ple poles for X = 0,a;, 2a;, . . . , these poles are precisely 
the eigenvalues of the uncoupled bosonic mode. Then 
the regular energy spectrum of the Rabi model in each 
invariant subspace 'H± with parity ±1 is given by the ze- 
roes of G±{x): for all values x^ with G±{x^) — 0, the 
n-th eigenenergy with parity ±1 reads = x^ — g^ jio. 
The functions G±{x) are plotted in Fig H] for A = 0.4 
and Lo = g ~ \ between x — ~\ and x — h. Their ze- 
roes determine the first six eigenenergies with parity -1-1 
and the first five (including the groundstate) with parity 
— 1. For special values of model parameters .g, A, there 
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FIG. 1: G+(x) (red) and G^x) (blue) 
for g = 0.7, a; = 1 and A = 0.4 



in the interval! [—1,5] 



are eigenvalues which do not correspond to zeroes of ([3|) , 
these are the exceptional ones. All exceptional eigenval- 
ues have the form E'^ = nw — g"^ /lu, that is, they lie on 
one of the so-called baselines [21[ and coincide with some 
point in the spectrum of the limiting case A = 0, which 
corresponds in spin-boson language to zero hybridization. 
All exceptional eigenvalues are doubly degenerate with 
respect to parity. The necessary and sufficient condition 



for an exceptional eigenvalue to lie on the n-th baseline 
reads, 



Kn{n) = 0, 



(6) 



which furnishes a condition on the model parameters g 
and |A|. As all the Kn{n) are independent polynomials 
m A^, there can exist at most two exceptional eigenvalues 
for given g,A. These exceptional solutions to (P) have 
been known for a long time and were first discovered 
by Judd [2l[. They occur when the pole of G±{x) at 
Xn — nuj is lifted because its numerator in Q vanishes 
which happens only if ([5]) is satisfied. 
The functional form of G±{x) reads. 



G±(x) 



Gl 



(x) 



(7) 



where Gj.{x) is entire in x. The position of the solutions 
to G±{x) = is dictated by the pole-structure of G±{x), 
which leads to the conjecture that the number of eigen- 
values in each interval [noj, {n + 1)lu] is restricted to be 0, 
1, or 2. Moreover, for large energies x ^ g. A, the entire 
part of G± (x) can be approximated by 



G°±(x)=(lT- 
' X 



exp 



V 2uJ 



(8) 



which is monotonous, always > and slowly varying on 
the scale set by to. This suggests that the relative po- 
sition of the solution Xj{n) in the vicinity of the n-th 
baseline is fixed by the sign of alone: If ft.^ > 
[hfi < 0) then Xj{n) lies in the interval [(n — l/2)w,na;] 
{[nu, (n + l/2)cj]). It would also entail that an interval 
[noj, (n -|- l)a;] with two roots of G±{x) = can only 
be adjacent to an interval with one or zero roots; in 
the same way, an empty interval can never be adjacent 
to another empty interval. These conjectures about the 
distribution of the large eigenvalues, which can be con- 
firmed numerically, lead to a fairly regular distribution of 
the energies and a violation of the Berry- Tabor criterion 
[l^ll^. Figl2] shows the lowest part of the spectrum for 
A = 0.4, UJ = I and g between and 1. One may prove 
that there are no level crossings within each parity sub- 
space, allowing the unique labeling of each state with 
a pair of two quantum numbers, \ip) = |no, ni): The par- 
ity quantum number no, which takes the values -1-1 and 
— 1, and ni — 0,1,2,..., which denotes the ni-th zero 
of Gno (x) . The exceptional solutions correspond to level 
crossings between T-L+ and "H-. This characterization 
of each eigenstate through two quantum numbers corre- 
sponding to the degrees of freedom of the system parallels 
the unique assignment of three quantum numbers n, Z, m 
to the eigenstates of the hydrogen atom, refiecting the 
quantization of radial and angular degrees of freedom, a 
hallmark of integrability. 

It seems therefore natural to call a quantum system in- 
tegrable when such an assignment can be made — inde- 
pendent of the explicit determination of conserved quan- 
tities or even action variables, which is only possible if 
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FIG. 2: Rabi spectrum for A = 0.4, = 1 and < p < 0.8 
in the spaces with positive (red) and negative (blue) parity. 
Within each space the states are labeled with ascending num- 
bers 0, 1,2,.. .. This labeling does not change with g because 
no lines intersect within spaces of fixed parity. However, level 
crossings of states with different parity occur. The spectral 
graph is composed of two intersecting "ladders" of level lines, 
each corresponding to one parity subspace. This labeling is 
used on the right side of the figure. On the left side the states 
with g — are labeled by the uncoupled degrees of freedom, 
i.e. in |±,n), +/— corresponds to the two-level system and 
n = 0, 1, 2, ... to the eigenstates of the bosonic mode. 



the system under consideration has an integrable classi- 
cal limit in the sense of Liouville [131 ■ Without making 
the assumption of a classical limit, our criterion reads. 

Criterion of quantum integrability: If each eigenstate of 
a quantum system with fi discrete and /2 continuous 
degrees of freedom can be uniquely labeled by /i + /2 = / 
quantum numbers {di, . . . , d/^ , ci, . . . c/j}, such that the 
dj can take on dmi{Hj) different values, where Hj is the 
state space of the j-th discrete degree of freedom and 
the Ck range from zero to infinity, then this system is 
quantum integrable. 

The criterion does not presuppose the existence of a fam- 
ily of commuting operators whose different spectra are 
associated with the {dj,Ck} but provides a condition on 
the spectral graph of the system, that is, the spectrum as 
function of a parameter, typically one of the coupling con- 
stants. Without such a deformation parameter (which 
must conserve integrability) the association of more than 
one quantum number to the levels of a non-degenerate 
spectrum is ill-defined and would be restricted either to 
models solvable via Bethe ansatz 25| or systems with 
integrable classical limit. As is well-known, the Berry- 
Tabor criterion [22| relies precisely on the existence of 
this limit, that is, it applies only to the correspondence 
limit of large quantum numbers and presumes the valid- 
ity of the semi-classical quantization scheme. The crite- 
rion proposed here is not restricted to the upper part of 
the spectrum or the existence of a classical limit. More- 
over, it can be given a phenomenological formulation us- 
ing the spectral graph: If the system with coupling g has 
/ > 1 discrete or continuous degrees of freedom, each 



eigenstate can be written as \ni, . . .rif, g). Upon varying 
the coupling g, each of these states defines an energy level 
E{ni, . . .nf;g) as function of g. The totality of states for 
all sets {rti, . . .Uf} forms an /-dimensional manifold of 
spectral lines, which will typically intersect each other 
when energies corresponding to different invariant sub- 
spaces become accidentally degenerate at special values 
of the coupling. In a fully nondegenerate spectrum, a 
single integer quantum number which assigns 0,1,2,... 
to the eigenstates with ascending energies would suffice 
for a unique description of each state — but whenever a 
level crossing occurs, this is no longer possible. A sec- 
ond quantum number is necessary to discern the ener- 
getically degenerate states. The crossing is accidental in 
the sense that no new symmetry appearing at the given 
value of the coupling parameter is responsible for the 
degeneracy, which will concern typically only two states 
in the spectrum. It is merely the fact that for these 
states, belonging to dynamically decoupled subspaces, 
the energy dependence on the coupling coincides. The 
second quantum number labels the invariant subspaces 
and can be used to discern the energetically degenerate 
states at the crossing point. According to the proposed 
criterion, Integrability is equivalent to the existence of 
/ numbers to classify energy levels uniquely, if the sys- 
tem has / continuous or discrete degrees of freedom. It 
should be emphasized that these "quantum numbers" are 
a more general concept than the radial and angular quan- 
tum numbers known from atomic physics — they need 
not correspond to physical quantities (actions) quantized 
in integer multiples of H and have nothing to do with 
the Bohr-Sommerfeld semi-classical quantization rules. 
They are only defined with respect to the unique de- 
scription of eigenstates; the integrable systems differ from 
the non-integrable ones because they allow for a "fine- 
grained" description through an /-dimensional vectorial 
label, whereas the latter have a one-dimensional label 
corresponding to energy as the only conserved quantity. 
In the Rabi model we have /i = /2 = 1 and degeneracies 
takes place between levels of states having different par- 
ity, whereas within the parity subspaces no level crossings 
occur. The spectral graph consists of two "ladders" of 
level lines |no, m) for m = 0, 1, 2, . . .. Each ladder corre- 
sponds to an invariant subspace of the Z2-symmetry char- 
acterized by no = ij the parity eigenvalue. Within each 
subspace the system corresponds to a single continuous 
degree of freedom and is therefore integrable by defini- 
tion. The global label (valid for all values of the coupling 
g) is two-dimensional as / = /i -I- /2 = 2; the Rabi model 
belongs therefore to the class of integrable systems. It 
may be, however, that the symmetry is even stronger 
than necessary to achieve integrability, analogous to the 
classically "superintegrable" systems [g^. The Jaynes- 
Cummings model is an example for this case. Here, the 
continuous [/(l)-symmetry leads to a further decomposi- 
tion of the subspaces with fixed parity into a direct sum 
of two-dimensional invariant spaces labeled by the un- 
bounded quantum number C — 0, 1, 2, . . . ; even (odd) 
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values of C correspond to odd (even) parity. The larger 
number of dynamically decoupled state spaces entails a 
second possibility to label the states uniquely: through 
C and a two- valued index tiq — ±. Now, all levels with 
different C may intersect, leading not to just two but 
infinitely many ladders labeled by C, each having two 
rungs, labeled by tiq. In Figl3]the four lowest levels of the 
Rabi model with positive parity (red lines) are compared 
with the corresponding levels of the Jaynes-Cummings 
model with C = 1,3,5 (black lines). The enlarged sym- 
metry of the latter leads to two level crossings which are 
not present in the Rabi model [27| . The appearencc of in- 
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FIG. 3: The spectrum of the Jaynes-Cummings model (black) 
compared with the Rabi model (red) for A = 0.4 and even 
parity. The state labeling of the former has the form li/)) = 
\no,C) with no = ±. Two accidential crossings occur between 
levels with C = 5 and C = 3at(;ft;0.5 and between C = 1 
and C = 3 at g ^ 0.73. 

tersecting ladders in the spectral graph can be detected 
without knowledge of the exact solution or the correct 
assignment of quantum numbers to the different levels, 
the only condition being a sufHcient numerical resolution 
to discern degeneracies from narrow avoided crossings. 
This is a phenomenological virtue of the proposed cri- 
terion which could be used in computer experiments to 
test whether a given numerically solvable system pos- 
sesses a hidden integrable structure. Although a large 
number of level crossings as a function of model param- 
eters gives a strong hint to integrability, it is difficult to 
make the argument quantitative, because the number of 
intersecting ladders could become infinite already for two 
continuous degrees of freedom. For systems within the 
present class of one continuous and one discrete degree of 
freedom, however, its application is fairly obvious. The 
proposed criterion is sufficient but not necessary: in sys- 
tems with factorized scattering matrices, a unique asso- 
ciation of eigenstates with quantum numbers is possible 
which treats the discrete degrees of freedom differently 
from the present scheme, via the so-called nested Bethe 
ansatz [2^. But also this labeling leads to intersecting 
ladders as function of the coupling parameter. 

On the other hand, the absence of any level crossings 
in the spectral graph is sufficient for non-integrability if 
the total number of degrees of freedom (continuous and 



discrete) exceeds one: it means that the states can be 
classified only by energy, the single conserved quantity 
always present in hamiltonian systems, and no invariant 
subspaces exists. In the context of the Rabi model this 
case can be realized by breaking the Z2-symmetry. A 
possible generalization of ([Ij reads 

= Loa) a + gax{a + a)) + ea^ + Actz. (9) 

The term eux breaks the parity symmetry which couples 
the bosonic mode and the two-level system. Physically it 
corresponds to a spontaneous transition of the two-level 
system which is not driven by the radiation field. The 
state space does not separate into two subspaces and in- 
deed the spectral graph exhibits no level crossings at all if 
e is not a multiple of w/2 This is shown in FigU) In 
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FIG. 4: The spectrum of the generalized Rabi model with 
broken Z2-symmetry and A = 0.7, e = 0.2. The two-fold 
labeling of the states corresponds to the two integrable limits 
g = on the left and <; — >■ oo on the right. For finite g neither 
labeling classifies the states properly. 

this situation, the eigenstates can be uniquely numbered 
as belonging to the n-th energy level in ascending order, 
\ip) — \n). We have only one quantum number, energy, 
corresponding to the sole conserved quantity. Because 
the number of degrees of freedom nevertheless exceeds 
one, this model must be considered non-integrable. In- 
terestingly, ^ is still exactly solvable, although it does 
not possess any symmetry. Define the functions 

oo 

R^{x) = Y,Kt{x)[^Y (10) 

n=0 

and 

RH^) = t^^[^T- (11) 

n=0 

The Kf^ [x) are again recursively defined, 

nKt = /,ti(x)i^±_i - Kt2, (12) 
with the initial condition K'^ = l,K^{x) = f^{x) and 

fti^)-- + ^(nu^-x±e+ (13) 

oj Zg \ X — nu! ± e y 
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The n-th eigenvalue En of ([9]) is given by the n-th zero 
x„ of 



Ge(x) = A^R+{x)R-{x) - R+{x)R-'{x) 



(14) 



through En = x„ — g^/w. The fact that can be diago- 
nahzed analytically although not even a discrete symme- 
try is present signifies that integrability and solvability 
are not equivalent in the realm of quantum physics. In 
contrast to classical mechanics, non-integrable quantum 
systems with exact solutions exist. 



II. CONCLUSIONS 

We have seen that the Rabi model is integrable con- 
trary to common belief. The exact spectrum is given 
by the zeroes of a well-defined transcendental function 
G±{x) which converges everywhere except at the set of 
points X G {nuj\n = 0, 1, 2, . . .}, where it has simple poles. 
The previously found exact eigenvalues for special values 
of the parameters g, A are always two-fold degenerate 
and correspond to the merging of a zero with a pole of 
G±{x) at a: = nuj. The eigenfunctions can be given ex- 
plicitely as well (see appendix) . The solution of the Rabi 
model appears in this way to be of comparable complex- 
ity with the well-known solution of the one-dimensional 
potential well, although it is a strongly interacting, fully 
quantum mechanical model. Its generalization with bro- 
ken Z2-symmetry is non-integrable according to our cri- 
terion on integrability, which is tailored to the present 
class of systems with more than one but less than two de- 
grees of freedom - nevertheless it is exactly solvable and 
describes probably the smallest non-integrable physical 
system. Whether exact solutions are possible for non- 
integrable quantum models with two or more continuous 
degrees of freedom is under current investigation. 



III. APPENDIX 

To derive the result ©-(O for G±{x) we use the repre- 
sentation of bosonic creation and anihilation operators in 
the Bargmann space of analytical functions in a complex 
variable z [231, 



a 



d_ 

dz 



(15) 



Then, after applying a Fulton- Gouterman transforma- 
tion and setting w = 1, the time-independent Schrodinger 
equation in the subspace with even parity reads 

^^^(^) + .9 + ^) ^(^) = ^VX^) - A^(-z), (16) 

which is a functional differential equation in z. The so- 
lution ^-!{z) is analytic in the whole complex plane if 
E belongs to the spectrum of Hb,- With the notation 



■0(z) = 4>i{z), ip{—z) = 4'2{z), one obtains a coupled 
system of first-order equations, 

[z + g)^M^) + {gz - E)Mz) + AMz) = 
[z - g)£Mz) - {gz + E)Mz) + AMz) = 0. 

With y ^ z + g, X = E + g'^, cj)i^2 = e-ff+f'^/^i^a, it 
follows, 



dy 



(18) 



{y~2g)—^P2 = ix-Ag^ + 2gy)4,2-A^Pi. (19) 
dy 

■02(2/) is expanded into a power series in y, 

oo 



n— — oo 



and from (fTS)) one obtains for ipi{y), 



A 



(21) 



n— — oo 



Ea. ([T9)) is then equivalent with the recurrence Q and 
([5]). For ^"1,2(2/) to be analytic in y (i.e. at the point 
z = —g) it is necessary that Kn{x) = for all n < 0. 
If the coefficient Kq is set to 1, this entails Ki{x) = 
fo{x), which fixes the initial conditions of However, 
the power series in y, though analytic at z = —g, has 
the finite radius of convergence R — 2g, which can be 
deduced from the asymptotic value l/(2g) of fn-i{x)/n 
for n — > 00. Therefore, ip2iz + g) will develop a branch- 
cut at z = (7, if the parameter x does not belong to the 
discrete subset of M which determines the spectrum of 
Hji in ?^_|_. But this condition on x follows easily from 
the fact that we have two representations for ip{z) in 7^+, 
one constructed via (j)2{z) and the other with 4'i{z), 

00 

V'(z) = 02(-z) ^ eS-^i^„(a;)(-z + .g)" (22) 

V'(z) = 0i(z) = e-s- Vi^„(x)A^i±^. (23) 
^ — ' x — n 

The expansion in powers of z is analytic at z = 5 in (j22p 
and analytic at z = —g in (j23p . Only if both expansions 
are analytic at z = ±g they may coincide everywhere and 
represent the same ?/'(^) for all z g C. Therefore 

G+{x;z) = <j)2i-z)~(t>i{z) = for aU zgC (24) 

if and only if a; = E + g'^ corresponds to a point in 
the spectrum of Hji. Because x is the only variable 
in G+{x;z) besides z, it suffices to solve G+{x;z) — 
for some arbitrarily chosen z. However, the function 
G+{x]z) is only well-defined in x via its expansion in 
powers of z within the joint radius of convergence of ( 
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and ([23)). which restricts the absolute value of z to be 
less then g. If solves G+(a;;z) = for one such z, 
it will solve (|24|) as well and the eigenfunction ■0(z; E) 
for the energy E — xq ~ possesses the two alternative 
series expansions (|22I23I) which converge for x = xq in 
the whole z- plane and are uniquely determined by xo- It 
follows that all regular energy levels of H]i with positive 
parity are non-degenerate. Setting now z = 0, we obtain 
G+{x) — G+{x;0). This argument for carries over 
to 'H-, the subspace with negative parity, by replacing 
A with —A, which completes the derivation of ([3]). An 
exceptional solution occurs when a zero of G± (x) merges 
with a pole at x = nuj. As then (j)2{—z) ^ 4>i{z), this 
point corresponds to a two-dimensional representation of 
Z2 and the eigenvalue is two-fold degenerate. 

There have been early attempts [29 31] to obtain 
an equation for the eigenvalue E{x) without the Z2- 
symmetry underlying the functional equation (|16p . as 
(|17p arises after a rotation in spin space directly from 
dD). Then, of course, no functional relation between 
and 02 can be used. Schweber has found the recurrence 
(|4l5p for the coefficients Kn{x) in 



(25) 



ri=0 



and noted that the solution belongs only to the Hilbert 
space if a branchcut aX y = 2g is avoided by the correct 
choice of x. One may formulate the problem in terms of 
the two kinds of solutions to the three-term recurrence (|4]) 
with initial condition Kq = 1. The dominant solutions 
[33 [ all have the asymptotics 



lim 



1 



Kn-i{x) 2g 



1 



(26) 



with finite radius of convergence R — 2goi ((25)) whereas 
the minimal solution has VrnVn^ao Kn{x) / Kn-i{x) = 0, 
which entails an infinite radius of convergence of (I^S]) . 
as mandated by analyticity. This imposes a second con- 
straint on the initial condition for Q besides 



Ko = 1, K^{x) = /o(x). 



(27) 



i.e. on the allowed value of Ki, because the minimal 
solution is unique, whereas the dominant solutions are 
only determined modulo a multiple of the minimal one. 
This second constraint fixes x and therefore the spec- 
trum. The procedure is then the following: 

1. ) Choose an arbitrary value x G R. 

2. ) Compute numerically the minimal solution to ^ 
with Kq ^ 1, {iv:™"'(x), n = 1,2,.. .}. The minimal so- 
lution depends on the fn{x) with n > 1. 

3. ) Compare /C{"™(x) with /o(x). 

If both coincide, x will correspond to the eigenenergy 
E = x-g^ 

In [29ll30j it is claimed that a formally exact solution 
of the eigenvalue problem may be thus obtained in the 
following way: Define Vn = Kn/Kn-i for n = 1,2,... 



Then ([4]) with initial condition (l27l) is equivalent to the 
two-term nonlinear recurrence relation 



fn~l{x) 



1 



nVn-l 



n>2 



(28) 



with initial condition Vi = foi^), if none of the Kn van- 
ishes. But (l28l) can also be written as 



1 



/„(x)-(n + l)K+i 



(29) 



which determines Vn "from above" instead by the initial 
condition Vi — fo{x) together with (|28l) . One could write 
Vi as a continued fraction. 



Vi{x) = 



1 



hix) 



(30) 



Because the r.h.s. of (PH)) depends only on the set 
/i(x), f2{x), . . ., the formal equation 



/o(x) = l|/i(x)-2|/2(x)- 31/3(2;) 



(31) 



seems to provide the sought condition on x, a transcen- 
dental equation, whose infinitely many roots determine 
the spectrum of Hj^. This argument, however, is falla- 
cious, because (|3T|) is only a reformulation of the recur- 
rence relation with initial condition Vi{x) = fo{x) 
and therefore valid for all x, regardless, whether they lie 
in the spectrum or not. Indeed, one may rewrite pi|) as 
the finite continued fraction 

/o(.t) = l\h{x) - 21/2(2;) - ... (n - l)|/„-i(x) - nVnix), 

(32) 

where Vn{x) depends on /o(a;) via psp and amounts then 
simply to 



fo{x) = foix). 



(33) 



Only if the r.h.s of (jSTj) is approximated in one way or the 
other, it may give a non-trivial equation for x, whereas its 
exact value is tantamount to p3|) . In other words, the 
continued fraction in p0|3ip has no well-defined value 
without specifying the limiting behavior of the Vn for 
n 00. It will converge both for the dominant solutions 
with Vn ^ 1/2(7 as well as for the sought minimal solution 
with Vn 0, the difference being that in the latter case 
the resulting Vi will not depend strongly on the cut-off 
procedure, i.e. setting Vn = s with e w for some 
sufficiently large N, whereas the choice of N and e has a 
strong impact on Vi , if one cuts off the continued fraction 
with Vn = e + l/2g [s^]. This, of course, follows from 
the uniqueness of the minimal solution to Q . 

Therefore pip is equivalent to the numerical algorithm 
described above. Moreover, the cut-off at n — N, neces- 
sary to define its r.h.s., renders it equivalent with a diag- 
onalization of Hn in a finite-dimensional Hilbert space, 
whose dimension depends on N. 

Clearly, this problem with the representation of the 
spectral condition on x in terms of an infinite continued 



7 



fraction is a consequence of neglecting the Z2-symmetry work was supported by the DFG through TRR80. 

which relates (jjii^) ^'^'^ 4'2{z) in a nontrivial way and 

allows the definition of the function G±{x), whose zeroes 

determine the energy eigenvalues of the Rabi hamilto- 

nian. 
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